Lipidation of Naturally Occurring α-Helical Antimicrobial Peptides as a Promising Strategy for Drug Design

In this paper, we describe the chemical synthesis, preliminary evaluation of antimicrobial properties and mechanisms of action of a novel group of lipidated derivatives of three naturally occurring α-helical antimicrobial peptides, LL-I (VNWKKVLGKIIKVAK-NH2), LK6 (IKKILSKILLKKL-NH2), ATRA-1 (KRFKKFFKKLK-NH2). The obtained results showed that biological properties of the final compounds were defined both by the length of the fatty acid and by the structural and physico-chemical properties of the initial peptide. We consider C8–C12 length of the hydrocarbon chain as the optimal for antimicrobial activity improvement. However, the most active analogues exerted relatively high cytotoxicity toward keratinocytes, with the exception of the ATRA-1 derivatives, which had a higher selectivity for microbial cells. The ATRA-1 derivatives had relatively low cytotoxicity against healthy human keratinocytes but high cytotoxicity against human breast cancer cells. Taking into account that ATRA-1 analogues carry the highest positive net charge, it can be assumed that this feature contributes to cell selectivity. As expected, the studied lipopeptides showed a strong tendency to self-assembly into fibrils and/or elongated and spherical micelles, with the least cytotoxic ATRA-1 derivatives forming apparently smaller assemblies. The results of the study also confirmed that the bacterial cell membrane is the target for the studied compounds.


Introduction
According to the report of the World Health Organization (WHO), the uprising resistance of pathogenic bacterial strains to commonly used antimicrobial drugs is claimed as one of the main global problems threatening public health and the economy [1]. Taking into account that antibiotic resistance is a natural phenomenon, it seems highly plausible that every subsequent antibiotic approved for treatment can induce, sooner or later, the emergence of resistant strains [2]. However, the misuse and/or abuse of antibiotics as well as inadequate sanitation and lack of effective prevention of infections accelerate the spread of microorganisms and the emergence of resistance to antibiotics [1]. Pathogenic microorganisms evade the effects of antibiotics by one of four mechanisms: limitation of drug uptake, modifying drug target, drug inactivation and active drug efflux [3].
Antimicrobial peptides (AMPs), due to their important role in the innate immunity of all living organisms, could serve as highly promising solution to the problem of drug resistance of pathogenic microorganisms [4,5]. A wide range of antimicrobial activity of AMPs against bacteria, viruses, fungi, and protozoa [6] predefines their potential use in the of their antimicrobial activity [24]. Previous studies have proven that the attachment of fatty acids to the N-terminus of peptides increases their activity against multidrug-resistant pathogens. However, the unbalanced extension of the hydrocarbon chain results in stronger interactions inside the formed aggregates and in the loss of antimicrobial properties [21,24]. The elongation of the fatty acid hydrocarbon chain is also correlated with the enhanced cytotoxicity of synthesized compounds against healthy human cells [24].
In the present work, we have studied the effect of lipidation on the antimicrobial properties of three naturally occurring α-helical peptides [25][26][27][28]: LL-I (VNWKKVLGKIIKVAK-NH 2 ), LK6 (IKKILSKIKKLLK-NH 2 ), and ATRA-1 (KRFKKFFKKLK-NH 2 ) (Table S1, Supplementary Material). The selected peptides fulfilled the following criteria: short amino acid sequence (up to 15 amino acids), high content of L-lysine residues, positive net charge, helicity, and confirmed activity against a wide range of microorganisms. Within the scope of our studies, we analyzed cytotoxicity, conformation and self-assembly properties of the most promising lipopeptides. In addition, we studied the interactions of the selected lipopeptides with artificial membrane models and visualized their effect on S. aureus cells. The performed studies allowed for a thorough analysis of the effect of lipidation on the structure-activity relationship.

Antimicrobial Activity Assays
Minimum inhibitory concentration (MIC) values against the representative strains of Gram-positive and Gram-negative bacteria, as well as fungi Candida spps, are presented in Table 1. As can be seen, the lipopeptides were active against Gram-positive bacteria and fungi, while Gram-negative E. coli was the most resistant strain. The antimicrobial properties were strongly dependent on the type of the attached fatty acid. The best activity with MIC values of 3.13-12.5 µg/mL against Gram-positive bacteria was observed for derivatives with hydrocarbon chain that comprised C 8 -C 12 carbon atoms. The extension of the carbon chain length with the consequent increase in the lipopeptide hydrophobicity led to reduction of the activity against representative Gram-positive bacterial strains but improved eradication of fungi strains, which is in a good agreement with the previous studies [29,30]. Nevertheless, in the case of palmitoyl-modified LL-I and LK6 analogues, a marked drop in antifungal activity was noted, while C 16 -ATRA-1 still remained active with MIC not exceeding 3.13 µg/mL. This discrepancy probably resulted from the highest overall positive charge of the ATRA-1 analogues. Vieira and Carmona-Ribeiro [31] showed that the critical step determining the antifungal effect of cationic compounds is the reversion of the cell surface charge from negative to positive rather than cell lysis, which well correlates with our results. Moreover, the highest overall charge decreased the ability of the peptide to selfassembly or facilitated dissociation of the monomers from assemblies during membrane association. Taking into account that the fungal envelope is more resistant to penetration by AMPs in comparison to the bacterial one, lipopeptide oligomerization could have a negative impact on antifungal properties, as it may impair the penetration ability of the lipopeptide through the cell wall and thus prevent interaction with the cell membrane [30].
Detailed analysis of the relationship between the optical density at 600 nm and the concentration of lipopeptides showed a surprising decrease in the antimicrobial activity in the concentration range exceeding the determined MIC values ( Figure S1, Supplementary Material). In a typical relationship, the pathogen's number decreases in a concentrationdependent manner over a whole concentration range. However, in the case of the studied lipopeptides, the observed anomaly arose from the lipopeptides' self-assembly, which reduced their affinity to the lipid bilayer of microorganisms [32].

Evaluation of Cytotoxicity
The results of the cytotoxicity of the most active lipopeptides against the immortal human keratinocytes cell line (HaCaT), human breast cancer cells (MCF-7 line, ATCC) and human prostate cancer cells (PC3 line, ATCC) are shown in Figure S2 (Supplementary Material) and Table 2. As can be seen, the compounds C 8 -ATRA-1 and C 10 -ATRA-1 are characterized by the lowest cytotoxicity to human keratinocytes. An increase in the length of the hydrocarbon fatty acid chain by two methylene groups reduced the cytotoxicity ca. 1.5 times. Both ATRA-1 conjugates showed high cytotoxic activity against the MCF-7 cell line, making them a promising candidate in anti-cancer drug development. The C 10 -ATRA-1 seems to be the most attractive, as its selectivity indexes are 9.47 (MCF-7) and 3.92 (PC3). This means that an almost 9.5-times lower lipopeptide concentration is required to inhibit the biological processes of MCF-7 cells and an almost 4-times lower concentration in the case of PC3 compared to healthy human cells. A similar tendency is observed for C 8 -ATRA-1. Unfortunately, in the remaining cases, lipidation gave a high cytotoxicity to both healthy and neoplastic cell lines.

Characterization of Peptide Assemblies
To examine a size distribution (hydrodynamic diameter) of oligomeric species of the selected lipopeptides, C 8 -LL-I, C 8 -LK6 and C 10 -ATRA-1, in water and phosphate-buffered saline (PBS, pH 7.4), dynamic light scattering (DLS) was performed ( Figure S3, Supplementary Material). The measurements were carried out for various concentrations of the lipopeptides (0.5-2 and 0.25-2 mM for salt-free and PBS solutions, respectively). The DLS results demonstrated the existence of multiple populations of self-assembled oligomeric species for the studied lipopeptides within the applied concentration range. In general, the average particle size distributions did not alter significantly with dilution. However, a decrease in concentration promoted the narrowing size distribution of nanostructures in most cases, except for C 10 -ATRA-1 in PBS. In 0.5 mM salt-free aqueous solution, the highest scattering intensity was centered at 404 ± 37, 327 ± 102 and 284 ± 105 nm for C 8 -LL-I, C 8 -LK6 and C 10 -ATRA-1, respectively. Transfer of the lipopeptides to the PBS resulted in a drastic shift of the peaks with the highest scattering intensity to larger values of the hydrodynamic diameters, 3360 ± 165, 2596 ± 512 and 2808 ± 1193 nm for 0.25 mM C 8 -LL, C 8 -LK6 and C 10 -ATRA-1, respectively. For comparison, the hydrodynamic diameters of the described compounds in 2 mM solutions in the PBS were found at 3844 ± 1031, 3047 ± 750 and 1040 ± 343 nm, respectively. Despite the apparent monodisperse distribution observed for C 8 -LL-I and C 8 -LK6, the spread of the hydrodynamic diameters reflected in the standard deviations is rather wide, indicating the inhomogeneous and fluctuating distribution of peptide assembly sizes. In the PBS solution, an increase in ionic strength partially screened the positive charges on the peptide molecules, reducing the electrostatic repulsions between them, and consequently facilitating the self-assembly. However, such large hydrodynamic diameters might also result from lateral association of lipopeptide assemblies and mutual diffusion. Therefore, to visualize the assemblies, the lipopeptides were examined with transmission electron microscopy (TEM) after 24 h incubation in the PBS solution. The TEM provided evidence for the formation of clumped and isolated fibrils, blocks of cylindrical micelles and spherical micelles, as shown in Figure 1.
saline (PBS, pH 7.4), dynamic light scattering (DLS) was performed ( Figure S3, Supplementary Material). The measurements were carried out for various concentrations of the lipopeptides (0.5-2 and 0.25-2 mM for salt-free and PBS solutions, respectively). The DLS results demonstrated the existence of multiple populations of self-assembled oligomeric species for the studied lipopeptides within the applied concentration range. In general, the average particle size distributions did not alter significantly with dilution. However, a decrease in concentration promoted the narrowing size distribution of nanostructures in most cases, except for C10-ATRA-1 in PBS. In 0.5 mM salt-free aqueous solution, the highest scattering intensity was centered at 404 ± 37, 327 ± 102 and 284 ± 105 nm for C8-LL-I, C8-LK6 and C10-ATRA-1, respectively. Transfer of the lipopeptides to the PBS resulted in a drastic shift of the peaks with the highest scattering intensity to larger values of the hydrodynamic diameters, 3360 ± 165, 2596 ± 512 and 2808 ± 1193 nm for 0.25 mM C8-LL, C8-LK6 and C10-ATRA-1, respectively. For comparison, the hydrodynamic diameters of the described compounds in 2 mM solutions in the PBS were found at 3844 ± 1031, 3047 ± 750 and 1040 ± 343 nm, respectively. Despite the apparent monodisperse distribution observed for C8-LL-I and C8-LK6, the spread of the hydrodynamic diameters reflected in the standard deviations is rather wide, indicating the inhomogeneous and fluctuating distribution of peptide assembly sizes. In the PBS solution, an increase in ionic strength partially screened the positive charges on the peptide molecules, reducing the electrostatic repulsions between them, and consequently facilitating the self-assembly. However, such large hydrodynamic diameters might also result from lateral association of lipopeptide assemblies and mutual diffusion. Therefore, to visualize the assemblies, the lipopeptides were examined with transmission electron microscopy (TEM) after 24 h incubation in the PBS solution. The TEM provided evidence for the formation of clumped and isolated fibrils, blocks of cylindrical micelles and spherical micelles, as shown in Figure 1.

Conformational Analysis
The CD spectra were recorded for the parent peptides and the most active lipopeptides in water, PBS and liposomes, POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)]) ( Figure 2). The zwitterionic POPC liposomes imitated eukaryotic membranes, while POPG liposomes with their anionic character imitated negatively charged microbial membranes [33]. In water and PBS, all the tested compounds adopted an unordered conformation. The addition of POPG or POPC lipids induced the formation of an α-helical structure with characteristic minima at 208 and 222 nm in the CD curves, except for ATRA-1 in POPC. As can be seen, lipidation induced visible increases in the helical content in the presence of the zwitterionic POPC phospholipids. Moreover, in the case of the LL-I and LK6 derivatives, a clearly higher helicity was observed in the POPC (40-48%) vesicles than in the POPG ones (19-27%). In turn, in the case of ATRA-1 derivatives, the opposite tendency was noted (Table 3), which correlated well with their lowest cytotoxicity. Moreover, the θ 222 /θ 208 ratio ≥ 1.0 determined for the parent LK6 and ATRA-1 peptides, and C 8 -LK6, C 8 -ATRA-1 and C 10 -ATRA-1 derivatives in the POPG, indicated the existence of a coiled-coil motif, which probably resulted from the formation of stable peptides' assemblies, while the θ 222 /θ 208 ratio ≤ 0.88 in the POPC for all studied lipopeptides suggested the presence of isolated α-helices [34][35][36]. Additionally, the conformational preferences of the selected lipopeptides, C8-LK6, C8-LL-I and C10-ATRA-1, were determined in the presence of S. aureus with various bacterial cell count of the inoculum ( Figure 2). The selection of S. aureus cells was determined by the results of antimicrobial assays, which clearly showed that this bacterial strain is the primary target for the tested compounds. Unexpectedly, the CD spectra indicated that C8-LK6 was the only one that adopted the α-helical structure (17-22%) in the presence of bacterial cells (Table 3). In less concentrated bacterial suspensions, the peptide tended to adopt a coiled-coil structure, but as the number of bacterial cells increased, the formation of an isolated α-helix was observed. With the increase in the bacterial cells/peptide ratio, the surface area of interaction of the lipopeptide with the cell membrane increases, which may affect the forces of interaction inside the coiled-coil structures and lead to the disintegration of peptide assemblies.

Effect of Lipopeptides on the Phase Transition Temperature of DPPC and DPPG Lipids
The changes in the structure of artificial membrane models, DPPC (phospholipids, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine) and DPPG (1,2-dipalmitoyl-sn-glycero-3-  Additionally, the conformational preferences of the selected lipopeptides, C 8 -LK6, C 8 -LL-I and C 10 -ATRA-1, were determined in the presence of S. aureus with various bacterial cell count of the inoculum ( Figure 2). The selection of S. aureus cells was determined by the results of antimicrobial assays, which clearly showed that this bacterial strain is the primary target for the tested compounds. Unexpectedly, the CD spectra indicated that C 8 -LK6 was the only one that adopted the α-helical structure (17-22%) in the presence of bacterial cells (Table 3). In less concentrated bacterial suspensions, the peptide tended to adopt a coiled-coil structure, but as the number of bacterial cells increased, the formation of an isolated α-helix was observed. With the increase in the bacterial cells/peptide ratio, the surface area of interaction of the lipopeptide with the cell membrane increases, which may affect the forces of interaction inside the coiled-coil structures and lead to the disintegration of peptide assemblies.

Effect of Lipopeptides on the Phase Transition Temperature of DPPC and DPPG Lipids
The changes in the structure of artificial membrane models, DPPC (phospholipids, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine) and DPPG (1,2-dipalmitoyl-sn-glycero-3-[phospho-rac-(1-glycerol)]) induced by the lipopeptides were monitored on FTIR spectra by tracking the position of symmetric CH 2 stretching vibrations of the lipid methylene groups, v s (CH 2 ), as a function of temperature ( Figure 3). In all cases the lipids underwent a phase transition from an ordered gel phase to a liquid crystalline phase, which was characterized by un upward shift of v s (CH 2 ) bands by ca. 3 cm −1 . The inflection point in the v s (CH 2 ) vs. temperature corresponded to a main phase transition temperature, T m . In the case of neat DPPC and DPPG membrane models, T m was found at ca. 39.5 • C. Premixing of the DPPC and DPPG vesicles with lipopeptides induced T m changes ranging from −2 to 1.5 • C. However, with C 8 -LL-I in DPPC, a shoulder observed at 45 • C on the first derivative curve indicated an increase in the T m even by 5.5 • C. The appearance of more than one maximum on the first derivative curve as well as the broadening of the phase transition suggest differential distribution of the peptides in the lipid bilayer.
With negatively charged DPPG, a rise in the T m temperature is provoked by electrostatic interactions between peptide and lipid headgroups, which reduce the repulsion between adjacent lipids, leading to closer contact between their acyl chains and consequently to their higher ordering. In turn, an insertion of molecules into the membrane driven by hydrophobic interactions typically leads to a decrease in the main phase transition temperature [37]. In the case of the studied lipopeptides, the relatively minor changes in the T m indicated more or less compensation for the electrostatic and hydrophobic interactions. This was particularly evident in the case of ATRA-1 derivatives, which despite the highest overall positive charge, did not induce changes in the T m of the DPPG or induced changes that were negligible below the resolution of the experiment. In turn, when comparing C 8 -LL-I with C 10 -LL-I, it was observed that elongation of the hydrocarbon chain of the fatty acid shifted the T m of DPPG to a lower value. This phenomenon confirmed an increase in the contribution of the hydrophobic interaction in the binding process. Additional information was deduced from the absolute values of the CH 2 stretching vibrational frequencies.
As can be seen in Figure 3, the lipopeptides initiated a shift of the CH 2 stretching vibrations to lower wavenumbers over the entire temperature range compared to the neat DPPG. This phenomenon, together with an elevated T m , confirmed a higher order of the lipid acyl chains and/or enhanced interchain vibration coupling. In turn, with the DPPC complexes, the symmetric CH 2 stretching bands were shifted to higher wavenumbers in both the gel and liquid-crystalline phases, which is usually related to higher membrane fluidity and/or reduced interchain vibrational coupling [38]. Both phenomena reflect disturbances in the lateral interactions between phospholipid chains in the complexes either due to increases in the number of gauche conformers in the hydrocarbon chains or their looser packing. In addition, the lipopeptides induced phase separation upon binding to the DPPC, leading to the formation of the lipid domains with different melting temperatures. The increase in the frequency of the CH 2 stretching vibrations over the whole temperature range correlated well with the existence of the lipid domains with reduced T m . Both findings confirmed increased fluidity of the DPPC lipid bilayer. On the other hand, the appearance of the domains with the T m above that of the neat DPPC implied a stabilization of gel-phase, which may be due to the dehydration of the membrane interface at the peptide binding site [39,40]. Moreover, above the T m , the choline groups of the DPPC are folded inward toward the surface, exposing more of the negatively charged phosphatidyl groups [41]. Thus, the possibility of electrostatic interactions between the positively charged lipopeptides and negatively charged phosphatidyl groups is rising. [phospho-rac-(1-glycerol)]) induced by the lipopeptides were monitored on FTIR spectra by tracking the position of symmetric CH2 stretching vibrations of the lipid methylene groups, vs(CH2), as a function of temperature ( Figure 3). In all cases the lipids underwent a phase transition from an ordered gel phase to a liquid crystalline phase, which was characterized by un upward shift of vs(CH2) bands by ca. 3 cm −1 . The inflection point in the vs(CH2) vs. temperature corresponded to a main phase transition temperature, Tm. In the case of neat DPPC and DPPG membrane models, Tm was found at ca. 39.5 °C.
Premixing of the DPPC and DPPG vesicles with lipopeptides induced Tm changes ranging from -2 to 1.5 °C. However, with C8-LL-I in DPPC, a shoulder observed at 45 °C on the first derivative curve indicated an increase in the Tm even by 5.5 °C. The appearance of more than one maximum on the first derivative curve as well as the broadening of the phase transition suggest differential distribution of the peptides in the lipid bilayer. With negatively charged DPPG, a rise in the Tm temperature is provoked by electrostatic interactions between peptide and lipid headgroups, which reduce the repulsion between adjacent lipids, leading to closer contact between their acyl chains and consequently to their higher ordering. In turn, an insertion of molecules into the membrane driven by hydrophobic interactions typically leads to a decrease in the main

The Thermodynamic Parameters of Lipopeptide Binding to the Lipid Bilayer
The results of the isothermal titration calorimetry (ITC) studies showed no, or more likely, too weak interactions between the lipopeptides and POPC, undetectable by the ITC method. The exception was C 10 -LL-I ( Figure S4, Supplementary Material), which bound to the POPC vesicles with a negative enthalpy change (∆H < 0) and a binding constant of 3.8·10 5 M −1 (Table 4). A closer inspection of the experimental data confirmed a spontaneous (∆G < 0) and entropy-driven (|∆H| < |T∆S|) process. Entropic control of the binding process indicated that hydrophobic interactions play a predominant role in the association of the peptide with zwitterionic POPC liposomes. In turn, all the ITC experiments with POPG liposomes confirmed the interactions of the lipopeptides with anionic lipids but with different thermodynamic effects. Thus, the titration of POPG to C 8 -LK6 and both ATRA-1 derivatives began with endothermic heat pulses accompanying the initial injections, becoming exothermic after several steps. These biphasic isotherms may suggest that the first stage of interaction was dominated by a release of a large number of water molecules both from the lipid bilayer and the peptide surfaces, which resulted in a large entropy gain. In turn, the exothermic effects following subsequent titration steps may be attributed to the peptide-membrane electrostatic attraction and peptide folding. On the other hand, the peptide-membrane association can occur in parallel with other processes, such as a change in the lipid phase properties, increase in permeability of the membrane, micellization of the lipid bilayer or initial peptide aggregation [42]. Despite the biphasic isotherms, better fitting results were achieved using the "one set of site" binding model, which showed that the overall process was dominated by an unfavorable enthalpy change (∆H > 0) offset by a large, favorable entropy change (T∆S > 0). With C 8 -LK6, the ITC experiment repeated with a twice higher initial peptide concentration ( Figure S4, Supplementary Material) highlighted even more complexity of the peptidemembrane interactions. However, computing accurate thermodynamic parameters for multiple binding events of many different affinities was not possible. Different results were obtained for the LL-I analogues. Namely, the titration of POPG into the peptide solutions resulted in a heat release over the whole concentration range. Unexpectedly, the binding of C 8 -LL-I to POPG vesicles was enthalpy-driven (|∆H| > |T∆S|), which indicated a significant role of electrostatic and van der Waals interactions and hydrogen bonds in the membrane association. In turn, an unfavorable entropic contribution to the free energy of binding (T∆S = −3.10 ± 1.36 kcal/mol) may suggest an enhanced ordering of the peptide-POPG complex. Titration of the POPG into the C 10 -LL-I solution involved a two-step binding process. An analysis of the isotherm with the "two sets of sites" model confirmed that both processes were spontaneous (∆G < 0) and entropy-driven (|∆H| < |T∆S|), but the determined parameters were burdened with a large error.

Transmission Electron Microscopy Imaging
Transmission electron microscopy was employed to visualize the effect of the representative analogue, C 8 -LL-I, on the S. aureus cells (Figure 4). The TEM micrograph of the untreated S. aureus cell showed a typical morphology with a relatively thick and, without doubt, intact cell wall. In turn, incubation of the S. aureus cells with C 8 -LL-I induced thinning of the membrane and its disintegration, as evidenced by its irregular shape. The results unambiguously confirmed that the destruction of the bacterial membrane is the primary effect that underlies the biological activity of the tested lipopeptides.
representative analogue, C8-LL-I, on the S. aureus cells (Figure 4). The TEM micrograph of the untreated S. aureus cell showed a typical morphology with a relatively thick and, without doubt, intact cell wall. In turn, incubation of the S. aureus cells with C8-LL-I induced thinning of the membrane and its disintegration, as evidenced by its irregular shape. The results unambiguously confirmed that the destruction of the bacterial membrane is the primary effect that underlies the biological activity of the tested lipopeptides.

Self-Assembly via Molecular Dynamics Simulations
The α-helical peptides are promising self-assembling building blocks in a bottom-up nanomaterial design strategy. Inspiration comes from inter-helical interactions that play an important role in the proper folding of many proteins, especially those attached to or associated with membranes [43][44][45][46]. Going further, conjugation of the α-helical peptide with a fatty acid increases the hydrophobicity of the compound, enhancing a tendency to self-assemble in solution and in membranes. Therefore, to take a deeper look at the self-assembly at the molecular level, full-atom self-assembly simulations were performed for the representative lipopeptides, C 8 -LL-1, C 8 -LK6 and C 10 -ATRA-1 ( Figures S5-S7, Supplementary Material).
The lipopeptides started to self-assemble in the first steps of the simulations to form oligomers with a number of aggregations up to nine ( Figure 5). The number of oligomers stabilized after approximately 250 ns of MD simulations for C 8 -LL-I and C 8 -LK6. With C 10 -ATRA-1, large fluctuations in the number of oligomers were observed until the end of the 500 ns MD simulations.
Both C 8 -LL-I and C 8 -LK6 self-assembled into a fibrillar morphology, producing staggered packing of α-helices, which allowed for axial growth of coiled-coil oligomers ( Figure 6). The interstrand hydrophobic interactions included most of all the N-terminal fatty acid chain and amino acid side chains in the central region of the peptides ( Figure 5D). These interactions stabilized the α-helical structure of the C 8 -LL-I and C 8 -LK6 during the MD simulations. In turn, C 10 -ATRA-1, despite the presence of the Phe-Phe motif [47] in the sequence, which is considered sufficient to promote self-assembly, showed the lowest self-assembly tendency among the tested lipopeptides. On the other hand, the ATRA-1 derivative had the highest overall positive charge among the tested peptides, which is not conducive to self-assembly. Therefore, for C 10 -ATRA-1, a large fraction of molecules in the monomeric state was reported after the 500 ns MD simulations, which promoted peptide unfolding and about a two-fold decrease in the overall helical fraction compared to the initial configuration. Among the oligomers, the most popular were those composed of 2-4 units. The trimers and tetramers adopted a globular shape rather than a fibrous structure with major hydrophobic interactions occurring between the decanoyl chains of adjacent molecules ( Figure 5D). the representative lipopeptides, C8-LL-1, C8-LK6 and C10-ATRA-1 ( Figures S5-S7, Supplementary Material).
The lipopeptides started to self-assemble in the first steps of the simulations to form oligomers with a number of aggregations up to nine ( Figure 5). The number of oligomers stabilized after approximately 250 ns of MD simulations for C8-LL-I and C8-LK6. With C10-ATRA-1, large fluctuations in the number of oligomers were observed until the end of the 500 ns MD simulations. Both C8-LL-I and C8-LK6 self-assembled into a fibrillar morphology, producing staggered packing of α-helices, which allowed for axial growth of coiled-coil oligomers ( Figure 6). The interstrand hydrophobic interactions included most of all the N-terminal fatty acid chain and amino acid side chains in the central region of the peptides (Figure  [47] in the sequence, which is considered sufficient to promote self-assembly, showed the lowest self-assembly tendency among the tested lipopeptides. On the other hand, the ATRA-1 derivative had the highest overall positive charge among the tested peptides, which is not conducive to self-assembly. Therefore, for C10-ATRA-1, a large fraction of molecules in the monomeric state was reported after the 500 ns MD simulations, which promoted peptide unfolding and about a two-fold decrease in the overall helical fraction compared to the initial configuration. Among the oligomers, the most popular were those composed of 2-4 units. The trimers and tetramers adopted a globular shape rather than a fibrous structure with major hydrophobic interactions occurring between the decanoyl chains of adjacent molecules ( Figure 5D).

Peptide-Membrane via Coarse-Grained Molecular Dynamics Simulations
Coarse-grained molecular dynamics simulations were used to investigate the interactions of the representative lipopeptides, C8-LL-1, C8-LK6 and C10-ATRA-1, with the S. aureus bilayer model. To validate the CG MD simulations, the results were compared to the intact lipid bilayer (Figures S8-S11, Supplementary Material). In the first steps of the CG MD simulations, the lipopeptide monomers were attracted to the lipid bilayer or self-assembled into oligomers and/or micelles, which finally approached the membrane

Peptide-Membrane via Coarse-Grained Molecular Dynamics Simulations
Coarse-grained molecular dynamics simulations were used to investigate the interactions of the representative lipopeptides, C 8 -LL-1, C 8 -LK6 and C 10 -ATRA-1, with the S. aureus bilayer model. To validate the CG MD simulations, the results were compared to the intact lipid bilayer (Figures S8-S11, Supplementary Material). In the first steps of the CG MD simulations, the lipopeptide monomers were attracted to the lipid bilayer or self-assembled into oligomers and/or micelles, which finally approached the membrane surface anyway. A detailed analysis of the results showed an expected drop in the number of clusters at the beginning of the CG MD simulations due to the lipopeptide self-assembly ( Figure 7A). However, for the C 8 -LK6 and C 10 -ATRA-1 systems, the number of clusters started to increase again, first relatively quickly (up to 600 ns) and then slowly, reflecting their insertion and gradual dispersion of the monomers in the outer leaflet of the membrane. Finally, all the C 8 -LK6 molecules and about 80% of the C 10 -ATRA-1 molecules penetrated into the lipid bilayer. Incorporation of the peptide's micelle into the lipid bilayer began with its gradual opening from the membrane-bound site, which was accompanied by a distinct bulge in both leaflets of the membrane ( Figure 7B). This phenomenon excluded polar residues from the contact region and reduced the electrostatic energy barrier for mixing hydrophobic tails of lipopeptides and membrane lipids. After complete insertion, the peptide molecules rapidly dispersed throughout the membrane, and the membrane itself flattened out. Finally, only negligible curvature was observed. By contrast, the C 8 -LL-I worm-like micelle, created quickly at the beginning, remained bound to the membrane surface until the end of the 10 µs CG MD simulations. Surface-bound micelles induced an evident segregation of membrane lipids, leading to a distinct enrichment of the micelle binding site in cardiolipin (CDL2) (Figures S9-S11, Supplementary Material). Following this, the increase in CDL2 density promoted the negative membrane curvature, a phenomenon well seen in the system with C8-LL-I ( Figure 8A). This also translated into changes of area per lipid (APL) and thickness of the membrane at the binding site ( Figure 8B,C and Figure S9C Supplementary Material). A deep analysis even showed a 2.4-fold increase in APLs compared to the neat membrane and various changes in the membrane thickness ( Figures S8 and S9, Supplementary Material). The center of the membrane space directly beneath the micelle was characterized by reduced thickness, and its periphery by increased thickness, with a difference of up to 13 Å. The perturbated bilayer structure was also evident from the twodimensional lipid density map in the upper leaflet of the membrane. As can be seen, the peptide disrupted the membrane integrity (Figure 8). Surface-bound micelles induced an evident segregation of membrane lipids, leading to a distinct enrichment of the micelle binding site in cardiolipin (CDL2) (Figures S9-S11, Supplementary Material). Following this, the increase in CDL2 density promoted the negative membrane curvature, a phenomenon well seen in the system with C 8 -LL-I ( Figure 8A). This also translated into changes of area per lipid (APL) and thickness of the membrane at the binding site ( Figure 8B,C and Figure S9C Supplementary Material). A deep analysis even showed a 2.4-fold increase in APLs compared to the neat membrane and various changes in the membrane thickness ( Figures S8 and S9, Supplementary Material). The center of the membrane space directly beneath the micelle was characterized by reduced thickness, and its periphery by increased thickness, with a difference of up to 13 Å. The perturbated bilayer structure was also evident from the two-dimensional lipid density map in the upper leaflet of the membrane. As can be seen, the peptide disrupted the membrane integrity ( Figure 8). Regardless of whether the positively charged lipopeptides were fully inserted into the membrane or only bound to it in an aggregated form, they affected the electrostatic state of the outer membrane-water interface. Moreover, they contributed to the displacement of sodium ions and changes in their concentration in the extracellular space. All this should lead to disruption of the electrostatic potential across the membrane ( Figures S9-S11, Supplementary Material), and consequently to the disturbance of the behavior and function of bacterial cells.

Peptide Synthesis
All the peptides were synthesized on a Tenta Gel S RAM resin (0.22 mmol/g) according to previously described methodology [42]. After deprotection of the N-terminal part of the peptide chain with 20% piperidine in N,N-dimethylformamide (DMF), the modifications were performed using fatty acids with various carbon chain lengths (C 8 -C 16 ), (1-cyano-2-ethoxy-2-oxoethylidenaminooxy)dimethylamino-morpholino-carbenium hexafluorophosphate (COMU) and N,N-diisopropylethylamine (DIPEA) of molar ratio peptide:fatty acid:COMU 1:3:3 (pH 8-9). Coupling of fatty acid was carried out for 2 × 15 min. Cleavage of the final peptides from the resins was performed using a mixture of trifluoroacetic acid (TFA), triisopropylsilane (TIPS) and water (95:2.5:2.5 v/v/v) for 2 h at room temperature (Scheme 1). Then, the resins were drained, and the obtained solutions were concentrated. The crude peptides were precipitated by the addition of diethyl ether and centrifuged; the films were then dried in a stream of nitrogen gas, dissolved in water and lyophilized.

Peptide Synthesis
All the peptides were synthesized on a Tenta Gel S RAM resin (0.22 mmol/g) according to previously described methodology [42]. After deprotection of the N-terminal part of the peptide chain with 20% piperidine in N,N-dimethylformamide (DMF), the modifications were performed using fatty acids with various carbon chain lengths (C8-C16), (1-cyano-2-ethoxy-2-oxoethylidenaminooxy)dimethylamino-morpholino-carbenium hexafluorophosphate (COMU) and N,N-diisopropylethylamine (DIPEA) of molar ratio peptide:fatty acid:COMU 1:3:3 (pH 8-9). Coupling of fatty acid was carried out for 2 × 15 min. Cleavage of the final peptides from the resins was performed using a mixture of trifluoroacetic acid (TFA), triisopropylsilane (TIPS) and water (95:2.5:2.5 v/v/v) for 2 h at room temperature (Scheme 1). Then, the resins were drained, and the obtained solutions were concentrated. The crude peptides were precipitated by the addition of diethyl ether and centrifuged; the films were then dried in a stream of nitrogen gas, dissolved in water and lyophilized. Scheme 1. Solid-phase peptide synthesis.
The synthesized analogues were purified by the Waters RP-HPLC preparative system with Jupiter 4 µm Proteo column, 90 Å, 250 × 10 mm (Phenomenex, Torrance, CA, USA). The linear gradient (B in A) from 10 to 100% within 60 min (A: 0.1% solution of TFA in water, B: 80% solution of acetonitrile in A) and with a flow rate of 5 mL/min. The purity of the lipopeptides (at least 95%) was performed by the PR-HPLC analytic system with Jupiter 4 µm Proteo column, 90 Å, 250 × 4.60 mm (Phenomenex, Torrance, CA, USA). The linear gradients were from 1 to 80% B, 15 to 90% B and 50 to 100% B within 30 min and with a flow rate of 1 mL/min. The identity of final products was carried out by recording the molecular mass by mass spectrometry (MALDI-TOF, Biflex III, Bruker, Karlsruhe, Germany). Analytical data of all the lipopeptides are presented in Table S2, Supplementary Materials.

Antimicrobial Activity Assays
The evaluation of the antimicrobial activity of the synthesized lipopeptides was performed against the following reference bacterial and fungal strains: The synthesized analogues were purified by the Waters RP-HPLC preparative system with Jupiter 4 µm Proteo column, 90 Å, 250 × 10 mm (Phenomenex, Torrance, CA, USA). The linear gradient (B in A) from 10 to 100% within 60 min (A: 0.1% solution of TFA in water, B: 80% solution of acetonitrile in A) and with a flow rate of 5 mL/min. The purity of the lipopeptides (at least 95%) was performed by the PR-HPLC analytic system with Jupiter 4 µm Proteo column, 90 Å, 250 × 4.60 mm (Phenomenex, Torrance, CA, USA). The linear gradients were from 1 to 80% B, 15 to 90% B and 50 to 100% B within 30 min and with a flow rate of 1 mL/min. The identity of final products was carried out by recording the molecular mass by mass spectrometry (MALDI-TOF, Biflex III, Bruker, Karlsruhe, Germany). Analytical data of all the lipopeptides are presented in Table S2, Supplementary Materials.

Antimicrobial Activity Assays
The evaluation of the antimicrobial activity of the synthesized lipopeptides was performed against the following reference bacterial and fungal strains:  [48] and yeasts [49].
According to the procedures, Mueller-Hinton broth (MHB) was used as the working medium for the analysis of the bacterial susceptibility, while RPMI-1640 (with L-glutamine, without sodium bicarbonate) growth medium buffered with MOPS (morpholinepropanesulfonic acid) and supplemented with 2% glucose was used for the experiments with the fungi.
The minimum inhibitory concentration values (MIC) were determined by exposing the bacteria or fungi cultures with an initial inoculum of 5 × 10 5 and 5 × 10 3 CFU/mL, respectively, to the increasing concentrations of the serially diluted tested lipopeptide applied in the range of 6.25-50 µg/mL. The experiments were carried out for 24 h at 37 and 35 • C for the bacteria and fungi, respectively. Then, the absorbance was measured at 600 nm for each well of the 96-well plate using a microplate reader (Sunrise TM, Tecan Trading AG, Männedorf, Switzerland). Based on the obtained results, the degree of inhibition of the growth of microorganisms (expressed in %) in the presence of the tested compounds range was determined with respect to the control probes. Inoculated wells without peptide were treated as a positive control (100% growth), and the wells containing working medium only were defined as a negative control (0% growth). All measurements were repeated in triplicate.

MTT Assays
MTT assay was used to determine the cytotoxicity of the studied compounds toward breast cancer cells (MCF-7 line, ATCC), prostate cancer cells (PC3 line, ATCC) and human keratinocytes (HaCaT line, ATCC). The MCF-7 cell line was cultured in a RPMI 1640 medium (Gibco, ThermoFisher Scientific, Waltham, MA, USA), the PC3 line in F-12K medium (Gibco, ThermoFisher Scientific, Waltham, MA, USA) and the HaCaT line in the high-glucose DMEM medium (Gibco, ThermoFisher Scientific Waltham, MA, USA). All media were supplemented with 10% fetal bovine serum (Gibco, ThermoFisher Scientific, Waltham, MA, USA) and antibiotics (penicillin and streptomycin, Gibco, ThermoFisher Scientific, Waltham, MA, USA). Cells were seeded in 96-well plates at a density of 4000 per well and incubated at 37 • C in a 5% CO 2 atmosphere. Then, the cells were treated with tested derivatives in the concentration range of 0.1 to 200 µg/mL (or 0.1-100 µg/mL in the case of C 8 -LK6 due to its limited solubility) or with the same volume of the vehicle (H 2 O) alone as a control. The prepared plates were further incubated for 48 h. After this time, the aqueous MTT salt solution (3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide, VWR) in a concentration of 4 mg/mL was added (25 µL/well). After 3 h of incubation, the formazan product was dissolved in DMSO, and the absorbance was measured at 570 nm (the background absorbance was measured at 660 nm) using an EnSpire microplate reader (Perkin Elemer, Waltham, MA, USA). The control cell viability was assumed to be 100%. Three independent experiments were carried out in three repetitions. The statistical analysis of results and IC 50 values estimation were performed using GraphPad Prism 7 software (San Diego, CA, USA).

Liposomes Preparation
Liposomes made of POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] and POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) were obtained by hydration of the dry lipid film. For this purpose, the lipids were dissolved in a mixture of chloroform and methanol in a volume ratio of 4:1. The solvents were evaporated in the stream of nitrogen gas, and then, the samples were lyophilized to remove residual organic solvents. Afterward, the lipid film was hydrated with phosphate-buffered saline (PBS, pH 7.4) and shaken for 2 h at 60 • C. Subsequently, the obtained multilamellar vesicles (MLVs) suspension was subjected to several cycles of freezing in liquid nitrogen and thawing at 60 • C to reduce the size and layers of the liposomes. Large unilamellar vesicles (LUV) were obtained by extrusion of MLV liposomes through polycarbonate membrane filters with a pore diameter of 100 nm (Whatman International Ltd., Kent, UK) using a mini-extruder (Avanti Polar Lipids, Inc., Alabaster, AL, USA).
In the FTIR measurements, MLVs built with DPPC (1,2-dipalmitoyl-sn-glycero-3phosphocholine) and DPPG (1,2-dipalmitoyl-sn-glycero-3-[phospho-rac-(1-glycerol)]) were prepared by dissolving the lipids in a mixture of chloroform and methanol in a volume ratio of 2:1, dried and lyophilized. The obtained films were resuspended in a phosphate buffer (pH 7.4) and mixed for 2 h at 60 • C. After this, the samples were frozen and thawed 5 times. The prepared MLVs were blended with the lyophilized lipopeptides to obtain the desired peptide-to-lipid 1:10 molar ratio and incubated for 1 h at room temperature.

Circular Dichroism (CD)
CD experiments were carried out at 25 • C using a Jasco J-815 spectropolarimeter in the wavelength range 185-260 nm. The concentration of the lipopeptides was 0.15 mg/mL. The CD spectra were recorded in water, PBS (pH 7.4), 1.3 mM POPC and 1.3 mM POPG. The results are plotted as the mean molar ellipticity per one amino acid residue (MRME, deg cm 2 dmol −1 ) as a function of the wavelength (nm). Additionally, we performed a conformational analysis of the selected compounds in the presence of S. aureus cells. In this case, 0.15 mg/mL solutions of the tested peptides were prepared in the phosphate buffer (10 mM, pH 7.0) and inoculated with bacteria (5 × 10 6 CFU/mL), cultivated the same way as for the MIC experiment. The mixtures were then incubated for 2 h before recording the CD spectra.
The percentage of the α-helical structure was calculated based on the mean molar ellipticity at 222 nm using the formula: where θ 222 corresponds to the value of the mean molar ellipticity at a wavelength of 222 nm [50].

Isothermal Titration Calorimetry (ITC)
The thermodynamic parameters of the interactions of the lipopeptides with POPG and POPC liposomes were determined by isothermal titration calorimetry (ITC). Measurements were performed at 25 • C using an AutoITC microcalorimeter (MicroCal Inc., Northampton, MA, USA) with a 1.4491 mL sample and reference cells. The experiment consisted in the titration of a buffered peptide solution with a concentration of 0.05 and/or 0.07 mM with a buffered solution containing POPG or POPC LUVs with a lipid concentration of 1.3 and 1.8 mM (10.02 µL; 29 injections, 2 µL for the first injection only), as it was described previously [51]. The result of the titration was a thermogram showing a series of peaks corresponding to the titration steps. A titration curve was obtained after integration of the peak areas. The obtained isotherms were fitted to the "one set of sites" or "two sets of sites" models available in the ORIGIN 7.0 (MicroCal Inc., Northampton, MA, USA), which made it possible to determine the enthalpy (∆H), the binding constant (K ITC ) and the stoichiometry (N) of the interactions between the lipopeptide and the lipids. The following formula was used to calculate the entropy change (∆S) and the Gibbs free energy (∆G): where: factor 55.5 is the molar concentration of water, R is the gas constant (1.986 cal·mol −1 ·K −1 ), and T is the absolute temperature.

FTIR Measurements
Fourier transform infrared spectra were recorded on a IFS66 spectrometer (Bruker, Billerica, MA, USA) equipped with a DTGS detector. The samples were placed between two CaF 2 windows separated with a 50 µm Teflon spacer. The spectra were measured between 3000 and 2500 cm −1 at 2 cm −1 resolution as the average of 10 measurements each with 16 scans. The cuvette was thermostated to the appropriate temperature in the range of 27-46 • C by using a Julabo F33 thermostat (Seelbach, Germany), and the sample temperature was controlled with a CHY502 thermometer. The spectrum of the phosphate buffer solution at the appropriate temperature was each time subtracted from the final spectra. The wavenumber positions for symmetric and asymmetric CH 2 stretching modes of the lipid hydrocarbon chains were measured using a multiple Gaussian curve fitting procedure in OriginPro 2021 software. Finally, the frequency of the symmetrical CH 2 stretching modes was plotted against the temperature. For a more precise determination of the temperature of the main phase transition of the lipids, T m , the first derivatives of the relationships were calculated. The main phase transition temperature of the DPPC and DPPG was determined before and after binding of the lipopeptides.

Dynamic Light Scattering (DLS)
The hydrodynamic diameters of the lipopeptide nanoaggregates were examined by the Litesizer 500 (Anton-Paar GmbH, Graz, Austria) at 25 • C. The measurements were carried out at various lipopeptide concentrations ranging from 0.5 to 2 mM and from 0.25 to 2 mM, in water and PBS, respectively. Two-fold serial dilutions were used to obtain the sequential concentrations. Each sample was sonicated for 10 min followed by incubation for 1 h at 25 • C and transferred to the polystyrene cuvette. The detector was positioned at a side scatter position (90 • ). The analysis of the obtained results and report generation was performed with KalliopeTM software (Anton-Paar GmbH, Graz, Austria).

Transmission Electron Microscopy (TEM)
TEM was used to observe ultrastructural changes in S. aureus cells in the presence of C 8 -LL-I. Bacterial cells were incubated overnight in a Mueller-Hinton medium. The cells were centrifuged (8000 rpm, 15 min) and resuspended in fresh media to obtain a highdensity cell suspension (5 × 10 5 CFU/mL). C 8 -LL-I in a MIC concentration (25 µg/mL) were mixed with the cells and incubated overnight at 37 • C. Afterward, the cells were centrifuged (8000 rpm, 15 min) and washed three times with PBS. The bacterial cells were fixed with 2.5% glutaraldehyde buffered at pH 6.5 with 0.1 M sodium cacodylate (Polysciences, Warrington, PA, USA) for 6 h at room temperature, as described earlier [52]. Additional fixation was performed with a 1% osmium tetroxide solution (Polysciences, Warrington, PA, USA) for 2 h at 4 • C. The cells were centrifuged and resuspended in 0.1 M cacodylate buffer, centrifuged again and dehydrated with ethanol. The bacteria were embedded in Epon 812 resin (Sigma-Aldrich, Darmstadt, Germany) at room temperature. A Leica UC7 ultramicrotome (Wetzlar, Germany) was used to obtain ultra-thin sections (55 nm). As contrasting agents, lead citrate and uranyl acetate were used. Cells were visualized using a Tecnai Spirit BioTWIN (EFI, Hillsboro, OR, USA) microscope at 120 kV.
To obtain TEM micrographs of the formed nanostructures, peptide samples in the concentration of 2 mg/mL were incubated in PBS at 37 • C for 24 h. Then, 5 µL of an appropriate sample was applied on carbon-coated copper grids (300 mesh), and after absorption time, the sample was stained with 2% (v/v) aqueous uranyl acetate. The micrographs were obtained using the same device and conditions as described above.

The Self-Assembly Simulations
The all-atom molecular dynamic simulations were used to trace the self-assembly properties of the selected lipopeptides. The simulations were performed using the GPU/CUDAaccelerated implementation of PMEMD in AMBER 18 [53]. Non-standard fatty acid residues were modeled with XLEAP. The point charges were optimized by fitting them to the ab initio molecular electrostatic potential (6-31G* basis set, GAMESS 2013-ab initio molecular electronic structure program) [54] for two different conformations, followed by consecutive averaging of the charges over all conformations, as recommended by the RESP protocol [55]. The initial configuration was constructed by randomly placing 50 peptide molecules in a simulation box using PACKMOL [56]. The system was solvated and neutralized by adding chloride ions using the LEAP program available in the AMBER 18 package and minimized for 40,000 steps (steepest descent method) with weak harmonic constraints on the peptide molecules (10 kcal/mol Å) used for the first 10,000 steps. Afterward, the system was equilibrated using a three-step protocol: (1) 10 ns simulations with isotropic pressure coupling and a 0.5 fs time step, (2) 10 ns simulation with isotropic pressure coupling and a 1 fs time step, and (3) 500 ns simulations with isotropic pressure coupling and a 2 fs time step. Periodic boundary conditions were employed in all the simulations. Long-range electrostatic interactions were evaluated by the particle Mesh Ewald (PME) summation. A cut-off of 10 Å was used for van der Waals interactions. The SHAKE algorithm was used to constrain bonds involving hydrogen. The temperature was maintained using the Langevin coupling scheme with a friction coefficient of 1 ps −1 , whereas a Berendsen barostat maintained a reference pressure set to 1.0 bar. The analysis of the MD trajectories was performed by AMBER 18 [53] and GROMACS 2019.5 tools [57]. The inter-residue contact maps were created with Python-based scripts available as a public repository on GitHub (https://github.com/strodel-group/Oligomerization-State_and_Contact-Map, accessed on 14 November 2022). Visualizations were performed with VMD 1.9.4a38 [58].

CG MD Simulations of Spontaneous Peptide-Membrane Interactions
Molecular dynamics simulations were carried out using the GROMACS 2019.5 package [57] and the MARTINI coarse-grained force field [59,60]. The lipid bilayer was composed of 54% DPPG, 36% Lys-PG (lysyl-phosphatidylglycerol) and 5% CDL2 (cardiolipin) equally distributed between two membrane leaflets to mimic the membrane of S. aureus [61]. The system was built using the insane.py script available on the martini website (http://cgmartini.nl/, accessed on 1 July 2022). The system was solvated and neutralized by adding counterions using the GROMACS tools. The salt concentration for bulk solution was 100 mM NaCl. The system was energy-minimized and equilibrated with the stepwise lowered force constant of the harmonic restraints (from 200 to 10 kJ mol −1 nm −2 ) to fix the position of the headgroups of the membrane lipids during simulations. After equilibration, the system was subjected to isothermal-isobaric molecular dynamics (NTP) with a 10 fs time step, as suggested by Wigner et al. [62] The temperature was held at 323 K using the v-rescale temperature coupling. The pressure was treated semiisotropically at 1 bar using the Parinello-Rahman barostat with a coupling constant of τ p = 12.0 ps. The relative dielectric constant for explicit screening was 15. Coulomb interactions were treated using a reaction field and a cutoff of 11 Å. Periodic boundary conditions were applied in all directions. In the next step, the S. aureus lipid bilayer after 5 µs of CG MD simulations was used to build the systems with the C 8 -LL-I, C 8 -LK6 and C 10 -ATRA-1 lipopeptides. The lipopeptide molecules were placed randomly close to the outer leaflet of the membrane and neutralized by adding chloride ions. Details of the initial systems construction are given in Table S3 (Supplementary Materials). The peptide-membrane systems were minimized and equilibrated with the harmonic restraints (10 kJ mol −1 nm −2 ) to fix the position of the headgroups of the membrane lipids and peptides during the initial simulations. Afterward, the system was subjected to 10 µs CD MD simulations with the conditions mentioned above. However, to prevent peptide molecules drifting away from the outer membrane-water interface and binding to the inner leaflet due to the periodic boundary conditions, flat bottom harmonic restraints were applied between each peptide molecule and membrane.
The data analysis was performed with g_lomepro software [63] and standard tools of the GROMACS. Visualizations were created in VMD 1.9.4a38 [58].

Conclusions
The α-helical lipopeptides presented in this study were characterized by high antimicrobial activity against a diverse range of Gram-positive bacteria and fungi from Candida spps. Based on the determined MIC values, C 8 -LL-I, C 10 -LL-I, C 8 -LK6, C 8 -ATRA-1 and C 10 -ATRA-1 were selected as the most promising compounds, and the S. aureus and S. epidermidis bacterial cells were selected as their main target. The cytotoxicity assays showed that the C 8 -ATRA-1 and C 10 -ATRA-1 analogues had relatively low cytotoxicity against human healthy keratinocytes but high cytotoxicity against the MCF-7 cell line. This could make them useful in the design of novel drug candidates active against breast cancer.
The tested lipopeptides demonstrated a strong tendency to self-assembly into fibrils and/or elongated and spherical micelles, while the size of the aggregates correlated well with cytotoxicity to the healthy keratinocytes. Among the lipopeptides with the highest antimicrobial potency, the C 10 -ATRA-1 analogue self-assembled into the smallest aggregates and at the same time revealed the lowest cytotoxicity. The ATRA-1 analogue carries the highest positive charge, suggesting that this feature affects self-assembly and contributes to cell selectivity. The C 10 -ATRA-1 analogue were also found to adopt a coiled-coil structure in the presence of the artificial lipid membranes. It is believed that the peptides folded in this way can be characterized by lower cytotoxicity [64].
The results of our study also confirmed that the bacterial cell membrane is the target for the studied compounds. The lipopeptides were able to penetrate the bacterial membrane and affect its properties, and monomers seemed to be the main active species. However, in silico studies showed that even the peptide assemblies that remained on the membrane surface induced drastic changes in the membrane permeability.
In summary, the main goal of our research was to design a novel group of peptides that form supramolecular assemblies and have biological activity. We chose α-helical antimicrobial peptides as the initial structures that were further derivatized with various types of fatty acids. The combination of nature-mimicking approaches inspired by the coiled-coil motif and the conjugation of helical peptides with hydrophobic alkyl chains to obtain peptide amphiphiles generated a wide repertoire of supramolecular hierarchic self-assembling architectures. Due to the combination of the physicochemical, in vitro and in silico analyses, our study contributes to the understanding of the role of the self-assembly in biological activity of the AMPs, which is of key importance for the potential use of this class of compounds in medicine and the pharmaceutical industry. Moreover, in our opinion, the obtained results expand the basic knowledge about peptide self-assembly, which is still poorly explored.